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摘 要: 相对 轨道 状态 确定 是 小 卫星 编队 正常 工作 的 基础 和 重要 保障 . 针对 小 型 化 、 低 成 本 的 星 载 
GPS 接收 机 定 轨 精 度 较 低 的 情况 ,提出 一 种 利用 单纯 星 间 测 距 信 息 对 卫星 编队 相对 轨道 进行 修正 
的 方法 . 引入 星 间 测 距 信息 ,采用 扩展 卡尔 曼 滤 波 算法 提高 编队 的 相对 轨道 估计 精度 . 通过 仿真 验 
证 ,证 实 了 该 方案 的 有 效 性 . 
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Correction of the Relative Orbit States Based on Inter-Satellite 
Ranging for Micro-Satellites Formation 
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(7. Key Laboratory of Electronics amd Information Technology for Space Systems， 
NSSC, CAS, Beying 100190, China; 
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Abstract: Determination of the relative states is the primary and important guarantee for normal operation 
of micro-satellite formation. With the miniaturization and low-eost spaceborne GPS receiver, the accuracy 
of orbit determination is low. In order to solve this problem, a method of correcting the formation relative 
orbit is given with only inter-satellite ranging information used. Inter-satellite ranging information is intro— 
duced and the extended Kalman filter is used to enhance the estimation accuracy of the relative orbit. The 
numerical simulation results show the validity of this approach. 
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小 卫星 编队 是 由 若干 颗 物 理 上 相互 分 离 的 小 卫 
星 、 微 小 卫星 、 纳 卫星 或 皮 卫 星 协同 工作 , 通过 无 线 
组 网 技术 形成 一 颗 “虚拟 整体 卫星 ”" ,共同 完成 空 
间 任 务 . 小 卫星 编队 上 县 有 研发 成 本 低 、 运 行 风 险 小 、 
鲁 棒 性 强 等 优点 ,还 可 根据 任务 需要 灵活 改变 编队 


* 中 国 科 学 院 专项 课题 基金 资助 项 目 ( Y42115DC1S) . 
收 稿 日 期 :2015-0242 


的 构 型 和 指向 ”, 在 科学 探测 和 军事 等 领域 具有 广 
阔 的 应 用 前 景 . 

小 卫星 编队 的 星 间 相 对 轨道 状态 确定 是 队 形 保 
持 、 编 队 成 员 协 同 工 作 的 重要 保障 . 因此 , 星 间 相 对 
轨道 的 精确 确定 十 分 必要 . 目前 有 多 种 方法 用 于 星 
间 相 对 轨道 状态 确定 ,主要 包括 : 激光 测量 法 、 红 外 
测量 法 、 视 觉 相 对 导航 法 ”、GPS 测量 法 ”等 . 由 于 
激光 波束 窄 ,需要 额外 的 引导 系统 ,增加 了 系统 的 复 
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杂 程 度 , 不 利于 测量 系统 的 小 型 化 ; 红外 测量 法 测量 
精度 较 高 ,但 是 只 能 提供 测 角 信息 ,不 能 提供 测 距 信 
息 ,无 法 单独 使 用 ; 视觉 相对 导航 法 测量 精度 较 高 ， 
技术 成 熟 , 但 是 其 作用 距离 较 近 ,不 适合 中 远程 的 星 
间 相 对 导航 ; 基于 星 载 GPS 测量 系统 的 相对 导航 法 
被 认为 是 有 效 的 高 精度 相对 状态 测量 敏感 器 "， 
GPS 测量 法 的 事后 处 理 相 对 轨道 状态 估算 精度 较 
高 ,在 文献 日 中 ,采用 GPS 测量 系统 获得 的 相对 轨 
道 精 度 达 到 了 毫米 量 级 ,但 是 要 想 获 得 较 高 的 实时 
确定 精度 ,需要 复杂 的 数据 处 理 过 程 ,小 型 化 的 CPS 
接收 机 很 难 实现 高 精度 的 相对 轨道 确定 . 针对 以 上 
问题 ,本 文 提出 了 采用 星 间 测 距 信 息 对 星 间 初 始 相 
对 轨道 进行 修正 的 方法 来 提高 相对 轨道 确定 精度 . 
编队 中 的 每 颗 卫星 上 都 安装 有 星 载 GPS 接收 机 , 利 
用 GPS 接收 机 作为 测量 设备 ,可 以 估算 出 从 星 与 主 
星 的 绝对 位 置信 息 , 则 从 星 与 主星 的 相对 位 置信 息 
可 以 利用 绝对 位 置信 息 做 差分 得 到 . 但 是 此 时 ,从 星 
与 主星 的 相对 位 置 的 误差 限 是 绝对 位 置 误差 限 的 2 
们 ,误差 限 被 放大 ,这 是 任务 所 不 希望 看 到 的 . 同时 
GPS 接收 机 要 实现 小 型 化 、 低 成 本 ,这 就 使 得 其 定位 
精度 进一步 降低 . 因此 ,本文 在 GPS 测量 系统 初步 
确定 了 卫星 编队 的 相对 位 置信 息 后 ,利用 星 间 测 距 
信息 对 星 间 相 对 轨道 进行 修正 . 


1 动力 学 模型 


本 文 提 出 的 卫星 编队 构 型 为 一 个 主星 与 多 个 从 
星 构成 的 空间 特定 队 形 , 为 了 便于 研究 问题 ,只 考虑 
主星 和 其 中 的 一 个 从 星 的 相对 轨道 状态 确定 情况 ， 
其 他 从 星相 对 轨道 状态 可 采用 同样 方法 来 处 理 . 
1.1 坐标 系 定义 

描述 卫星 编队 相对 运动 时 需要 用 到 惯性 坐标 系 
和 相对 运动 坐标 系 , 为 了 方便 描述 其 相对 运动 和 建 
立 动力 学 模型 ,下 面 给 出 这 两 种 坐标 系 的 定义 . 

惯性 坐标 系 N: 坐标 系 原 点 Ov 与 地 球 地 心 重 
合 ,X\ 轴 指 向 春分 点 方向 ,Z\ 轴 垂 直 于 赤道 面 指向 
北极 ,六 轴 与 XX 轴 、Z\ 轴 满 足 右手 定 则 . 

相对 运动 坐标 系 H: 坐标 系 原 点 04 与 卫星 质 
心 重 合 ,X 轴 方 向 与 主星 指向 地 球 的 矢 径 一 致 , 罗 
轴 垂 直 于 Zi 轴 , 与 主星 运行 速度 方向 一 致 ,Zi 轴 与 
XX 轴 、Y 轴 满 足 右手 定 则 . 

在 描述 从 星相 对 位 置 时 ,以 主星 所 在 的 相对 运 
动 坐标 系 作为 参考 坐标 系 . 从 星 在 参考 坐标 系 中 的 


位 置 矢 量 可 表示 为 -= kk y 4d ,主星 与 从 星 间 的 
距离 为 


1 卫星 编队 参考 坐标 系 


Fig.1 Reference coordinate system of satellites 


formation reference coordinate System 


1.2 卫星 动力 学 方程 
环绕 地 球 轨道 运行 的 卫星 不 仅 受 到 地 球 球形 引 
力 的 作用 ,还 受到 其 他 摄 动力 的 影响 ,如 地 球 非 球形 
景 动 、 大 气 阻 力 摄 动 \ 日 月 引力 摄 动 、 潮 汐 摄 动 等 , 根 
据 牛 顿 力 学 原理 ,卫星 在 惯性 坐标 系 下 的 动力 学 方 
程 为 
R=-KR+a (1) 
Re 


式 中 ,R 为 卫星 到 地 心 距 的 大 小 ,R 为 卫星 地 心 距 矢 
量 ; 为 地 球 引 力 系数 ,a 为 卫星 受到 的 除 地 球 引 力 
外 其 他 所 有 摄 动 加 速度 之 和 . 
1.3 ”相对 运动 方程 
从 星 在 参考 坐标 系 中 的 位 置 为 (x,y,z) ,速度 为 
(*,y,z) . 主星 为 圆 轨道 ,从 星 与 主星 的 距离 远 小 于 
主星 的 轨道 半径 . 因此 ,卫星 编队 的 相对 运动 可 以 采 
用 Hil 方程 来 描述 
TX -2 -nx = a, 
上 + 2nx = a, (2) 
z+nz= 4a, 


式 中 ,n 为 主星 的 平均 角速度 


We (3) 


a,,;4,;4; 为 相对 摄 动 加 速度 . 由 于 两 载荷 距离 较 近 ， 
所 受 摄 动力 情况 基本 相同 ,因此 ,可 忽略 相对 摄 动 影 
响 吕 , 即 c. = a, = a, =0. 求解 式 (2) 可 得 相对 位 置 
,yz 的 解析 表达 式 , 对 相对 位 置 求 导 可 得 相对 速 
度 *,y;z 的 解析 表达 式 . 采样 时 间 为 >, 则 离散 化 的 
Hil 方程 为 
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Xi = DX, + W, (4) 
式 中 ,X; 为 时刻 的 从 星相 对 轨道 状态 , Wi 为 过 程 
误差 阵 ,@$ 为 状态 矩阵 . 


$=- pi |， 
pr Pp 
式 中 
4 -3cos( nr7) 0 0 
pi = a sin(n7) -~n7) 1 0 | 
0 0 cos(n7) 
sin( n7) 2(1 - cos(n7)) 0 
2(cos( nr) -1) 4sin(n7) _37 0 
0 0 sin( n7) 
3nsin( n7) 0 0 
pp = 区 cos(n7) -1) 0 0 | 
0 0 -msin( nz7) 


cos( PT) 2sin( n7) 0 
p,, = —2sin(n7) 4cos(n7) -3 0 
0 0 


cos( n7) 
2 测量 原理 


每 个 编队 成 员 上 都 安装 有 一 台 星 载 CPS 接收 
机 ,同时 ,在 卫星 上 安装 有 无 线 电 测 距 装置 , 可 以 测 
出 主 ` 从 星之 间 的 距离 . 首先 , 由 星 载 GPS 测量 系统 
确定 出 主星 和 从 星 在 惯性 坐标 系 下 的 绝对 轨道 位 
置 ,再 对 其 进行 差分 处 理 ,可 得 出 两 星 的 初始 相对 轨 
道 位 置 . 将 初始 相对 轨道 信息 作为 下 一 步 的 测量 值 ， 
加 入 星 间 测 距 信息 ,采用 扩展 卡尔 曼 滤 波 算法 对 相 
对 轨道 进行 进一步 修正 . 测量 原理 如 图 2 所 示 . 


GPS 届 星 


GPS 卫星 


主星 从 星 


图 2 卫星 编队 测量 原理 图 


Fig.2 Satellites formation measurement schematics 


2.1 GPS 测量 模型 
GPS 测量 系统 基于 载波 相位 的 测量 模型 为 
L(t) = p(t) +e(6T,(t) -6T'(t -7(t))) + 
AN + M(t) + el(i) (5) 
式 中 ,i 为 信号 接收 时 间 ,p:( 1) 为 接收 机 和 GPS 卫星 
的 几何 距离 ,ce 为 光速 ,57,(1) 为 接收 机 在 t 时 刻 的 
钟 差 ,r'( 0) 为 信号 传输 时 间 ,57 (1 一 7(i) ) 为 GPS 
卫星 在 信号 发 射 时 刻 1 -7(1) 的 钟 差 , 和 A 为 波长 ,N 
为 整 周 期 模糊 度 , M( 1) 为 其 他 系统 偏差 ,e( 1) 为 测 
量 噪声 . 
2.2 星 间 测 距 信息 
无 线 电 测 距 系统 可 以 测 出 主星 与 从 星 间 的 距 
离 , 星 间距 离 表 达 式 为 
r= Vx + + + v (6) 
式 中 , Vx +y +z 为 主 \ 从 星 的 几何 距离 ,vw 为 测量 
随机 误差 . 


3 滤波 算法 


主 ` 从 星 通过 GPS 接收 机 获得 测量 值 后 ,利用 
在 惯性 系 下 的 动力 学 方程 式 ( 1) 和 测量 方程 式 (5) ， 
采用 适当 的 滤波 算法 可 解 算出 绝对 轨道 信息 ,然后 
通过 作 差 处 理 即 可 获得 星 间 相 对 轨道 初步 估计 值 
XX =k yz x yy zz]', 再 加 入 星 间 
测 距 信 息 , 采 用 扩展 卡尔 曼 滤 波 算法 对 星 间 相 对 轨 
道 状 态 进 行 二 次 修正 . 

选取 从 星 在 参考 坐标 系 中 的 位 置 和 速度 为 状态 
量 X= k y z % 》 悦 状态 方程 采用 式 (4) .将 
基于 GPS 测量 系统 得 到 的 初始 相对 轨道 估计 值 和 
作为 观测 量 , 再 加 入 星 间 测 量 信息 7, 得 到 的 观测 量 为 

Pl)=E ad 7 2 了 (7) 
则 观测 方程 为 


pr = Xi + Vi (8) 


ZE 


xk + YE + 2 
可 简 记 为 :p; =h( 头 ,) + VV. 其 中 ,VV 为 测量 误差 阵 . 

相对 轨道 状态 确定 的 扩展 卡尔 曼 滤波 算法 计算 
过 程 如 下 : 
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s。 4S。 


入 


首先 给 定 大 时 刻 的 状态 量 x,, ,状态 的 一 步 预 
报 值 x ,1,1 为 


状态 估计 误差 协 方 差 阵 为 
Pi = 『 - KH,,] Pi 『 - KH,,] + 


wai 三 Px, (9) Ku Ri Ki (12) 
预报 误差 协 方差 阵 为 式 中 ， 
尼 = DP, + 人 (10) oh( ,1) 1 
计算 滤波 增益 K,，， 人 
K,,, = Puli ,Pn Hi . R,,] ” (11) 其 中 ， 
H = | 2 7 2 = 5 3 - 人 0] 
xX 十 Y 十 2 xX 十 Y +Z xX 十 Y 十 和 


状态 的 最 优 估计 值 x; ,11,1 为 


入 


Xi+Dk+l = bE + 下 bis ~ h(k+ 1 ,X11,) ] 
(14) 

式 中 ,x ,1,4 为 利用 时 刻 状态 量 的 一 步 预 测 值 ， 
,i441 为 +1 时 刻 状态 量 的 最 优 估计 值 ,y,,1 为 
+1 时 刻 的 测量 值 ,有 到,, ,为 测量 方程 的 关系 矩阵 ， 
Pi,; 为 时 刻 的 误差 协 方差 阵 ，Q,; 为 系统 噪声 序列 
方差 阵 ,R, ,为 量 测 噪声 序列 方差 阵 . 

给 定 相对 状态 初 值 和 和 初始 误差 协 方差 阵 
PP,, 即 可 采用 EKF 算法 对 编队 相对 轨道 状态 进行 
估计 . 


4 仿真 结果 及 分 析 


主星 与 从 星 的 轨道 高 度 相 同 , 均 为 7 400 km. 主 
星 的 离心 率 为 0°, 轨道 倾角 为 30°, 升 交点 赤 经 为 
10° ,近地点 幅 角 为 60°, 主 、 从 星 间 距离 约 为 1 km. 
卫星 编队 相对 轨道 状态 滤波 修正 过 程 的 原理 如 图 3 
所 示 . 

作为 卫星 编队 真 值 的 相对 轨道 信息 计算 时 考虑 
了 态 项 摄 动 影响 , 这样 更 接近 真实 运动 情况 ,而 在 
滤波 算法 中 采用 了 没有 考虑 任何 摄 动 影响 的 Hil 方 
程 . 星 间 无 线 电 测 距 精度 为 lcm. 首先 初步 给 出 误差 
较 大 的 从 星相 对 轨道 状态 . 然后 以 初始 相对 轨道 信 
息 作为 下 一 步 修 正 的 测量 值 ,加 入 星 间 测 距 信息 ,以 
式 (4) 作为 状态 方程 , 式 (8) 作为 观测 方程 ,采用 
EKF 滤波 算法 对 相对 轨道 信息 进行 二 次 滤波 修正 . 
取 初 始 估计 的 相对 轨道 状态 作为 滤波 状态 初 值 X,， 
初始 协 方差 阵 取 局 = diag{ 100,100,100,100,100， 
100} ;过程 误 差 方 差 阵 取 @ = diag{ 10“ ,10…， 
10“,10-"”,10-",10-"} ,测量 误差 方差 阵 取 R = 
diag{fl,1,1,10,10 10,10…}. 修正 前 和 修正 
后 的 从 星相 对 轨道 状态 误差 曲线 如 图 4 和 图 5 
所 示 . 


星 载 CPS 接 收 
机 


至 靳 GPS 接收 
机 


从 星 绝对 
轨道 确定 


图 3 卫星 编队 相对 状态 确定 原理 图 
Fig.3 Schematic diagram of the relative state 


determination for satellites formation 
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图 4 修正 前 后 从 星相 对 位 置 误差 对 比 图 


Fig.4 The relative position error contrast 
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图 5 修正 前 后 从 星相 对 速度 误差 对 比 


Fig.5 The relative velocity error contrast 


表 1 相对 轨道 状态 的 估计 误差 


Tab.1 Standard deviation of the relative orbit states 


oi/m o,/m o./m 
修正 前 9.8060 9.890 3 10.007 4 
修正 后 0.204 2 0.1947 0.062 5 
ol( m/s) ow ll m/s) Ol( m/s) 
修正 前 0.0101 0.010 1 0.009 8 
修正 后 “5.8857x10-4 5.293 1x10-+ 2.7278x10- 


从 仿真 结果 可 以 看 出 ,未 采用 星 间 测 距 信息 修 
正 前 ,从 星相 对 位 置 误 差 最 大 约 为 20 m, 相 对 速度 
误差 最 大 约 为 0.01 m/s, 加 入 星 间 测 距 信息 对 相对 
轨道 状态 修正 后 ,稳定 段 相 对 位 置 最 大 误差 小 于 
1 m, 相 对 速度 最 大 误差 小 于 1 x10”“ m/s. 表 1 给 
了 修正 前 后 各 轴 的 估计 精度 . 

从 表 1 中 可 以 看 出 加 入 星 间 测 距 信 息 修正 后 ， 
星 间 相 对 轨道 状态 的 估计 精度 明显 改善 . 该 方案 中 
测量 方程 中 共有 7 个 测量 值 ,下面 考虑 减少 测量 值 
的 数目 后 是 否 还 能 达到 修正 的 目的 . 在 测量 方程 中 
减少 3 个 相对 速度 ,测量 方程 变 成 4 维 ,其 他 条 件 不 


变 , 对 其 进行 仿真 . 下 面 给 出 测量 方程 维 数 减少 后 的 
仿真 结果 : 从 星相 对 位 置 估 计 精 度 r、= 0. 593 5 m， 
o，= 0.676 3 ma。 = 0. 130 2 m ,从 星相 对 速度 估 
计 精 度 o,。= 9.453 4 x 10“ my/syo，= 9.820 1 x 
10™ m/s,o,, = 3.9024 x10™ m/s. 

从 仿真 结果 中 可 以 看 出 , 减少 测量 方程 的 维 数 
后 , 星 间 测 距 信息 对 相对 轨道 状态 仍 有 修正 作用 ,但 
是 修正 后 的 精度 有 所 下 降 . 测量 方程 的 维 数 减少 可 
以 减少 一 定 的 计算 量 , 降 低 硬件 的 负担 ,具有 工程 实 
际 意 义 . 


5 结 论 


本 文 针对 小 型 化 的 星 载 GPS 系统 对 星 间 相 对 
轨道 状态 估计 精度 不 高 的 情况 ,提出 了 采用 星 间 测 
距 信息 对 其 进行 修正 的 方案 . 在 测 距 精 度 为 1 cm 的 
情况 下 , 进行 了 仿真 验证 ,仿真 结果 表明 采用 星 间 测 
距 信息 进行 修正 可 以 明显 地 提高 相对 轨道 状态 估计 
精度 . 同时 ,文中 还 讨论 了 测量 方程 维 数 减 少 的 情 
况 , 此 时 对 相对 轨道 状态 仍 有 修正 作用 ,但 是 修正 精 
度 有 所 下 降 . 
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综 上 所 述 : 在 双 曲 线 轨道 自主 导航 中 , 远离 近 点 
的 轨道 段 导 航 精 度 很 低 , 但 对 目标 天 体 特性 依赖 小 ， 
且 在 靠近 近 点 的 位 置 进行 轨道 确定 的 效果 理想 能 够 
满足 近 点 自主 制 动 任务 要 求 的 精度 . 在 椭圆 轨道 导 
航 中 ,合理 选择 计算 周期 对 精度 有 较 大 影响 , 滤波 的 
效果 不 明显 ,不 适用 于 近 圆 轨道 
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